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Abstract 



The matrix elements of the deformation potential of C70 are calculated by 
means of a simple, yet accurate solution of the electron-phonon coupling prob- 
lem in fullerenes, based on a parametrization of the ground state electronic 
density of the system in terms of hybridized orbitals. The value of the 

calculated dimensionless total electron-phonon coupling constant is A 0.1, 
an order of magnitude smaller than in Cqq, consistent with the lack of a su- 
perconducting phase transition in C70A3 fullerite, and in overall agreement 
with measurements of the broadening of Raman peaks in C70K4. We also 
calculate the photoemission cross section of C^q, which is found to display 
less structure than that associated with C^q, in overall agreement with the 
experimental findings. 

PACS: 74.70.Wz, 63.20.Kr, 33.60.-q, 61.48. +c 
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I. INTRODUCTION 



The discovery that some alkah-doped Cgo compounds are superconducting with 
a transition temperature as high as 40 K for the fuUeride CeoCss under pressure has 
triggered considerable interest in the study of fuUeride doped materials at low temperatures. 
Much effort has been concentrated in trying to understand whether such large values of Tc can 
be caused by the coupling of electrons to phonons. The electron-phonon coupling constants 
for Cqo have been estimated both from experiments and calculations §]. There exists some 
consensus on the fact that Ceo compounds are s— wave BCS-like superconductors driven by 
the coupling to selected intramolecular phonons . On the other hand, no calculation of this 
coupling exists in the literature for the case of C70. Since the associated alkali compounds 
display no superconductivity, with a smooth and continuous behavior of the resistivity down 
to about 1 K 0, it is an important issue to compare the electron-phonon coupling constant 
for C70 with those of Ceo- 

In the present paper we address this problem from a theoretical point of view, comput- 
ing the deformation potential of C70 and extracting the associated total electron-phonon 
coupling constant, A. It will be concluded that the resulting value of A is much smaller than 
in Ceo, essentially ruling out the possibility of a superconducting phase at low temperatures 
for doped fullerides built out of C70. 

This negative result, which can be understood in rather general, qualitative structural 
terms, as well discussed at the end of Sec. Ill, has two positive oucomes: first, it provides 
further evidence for the s— wave BCS-like mechanism at the basis of alkali-doped fullerides 
superconductivity. Second, it sheds light on the properties fuUerenes have to exhibit to be 
"good" building blocks of eventual organic superconductors, indicating the direction of the 
highly symmetric, small fullerenes (cf. also Refs. [|TO| , pJ]| ), a strategy which seems to be 



universally valid also for other compact systems, like the atomic nucleus ||12|| . 

In the production of fuUerene clusters by means of laser vaporization or arc discharge 
the yield of C70 can be optimized, by varying the experimental conditions up to a 

ratio of about 10% of C70 with respect to Ceo- This allows the production of sufficiently 
large amounts of C70 for the synthesis of C70 molecular crystals [0. Hence, the structural 
and electronic properties of C70 can be experimentally studied, and compared with those of 
Ceo- In analogy with the case of solid Ceo, C70 crystals can be doped with alkali metals, and 
the conductivity of the resulting compounds can be investigated. 

On the theoretical side, ab-initio methods based on the local density approximation 
(LDA) [|1^] have been applied to the study of the ground-state properties of C70. The 
ionic configuration, characterized by a D5/J geometry, is consistent with the NMR measure- 
ments |]T^ , and the electronic structure agrees reasonably well with UPS data. We adopt in 
this work the ground-state geometry of Ref. |jl6|]. The vibrational properties of C70 have been 
studied in [^, starting with the same geometry. A simplified model like the bond charge 
model (BCM) is nevertheless able to reproduce the experimental values for the optically 
active frequencies within 4% \T^\. Indeed, beyond the electronic and vibrational properties, 
ab-initio methods also allow to compute the coupling between these two degrees of freedom 
without the need of introducing any adjustable parameter. However, ab-initio methods in 
the case of large clusters, such as Ceo and C70, are computationally demanding and may not 
be particularly transparent. As an alternative, in order to extract the electron-phonon cou- 
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pling constants, it is possible to adopt a simplified numerical method which already proved 
to be satisfactory in the case of Cgo This method, albeit approximated, has been shown 
to lead to quite satisfactory results, also because it is devised to take into account accurate 
phonon eigenvectors, which are known to influence very strongly the matrix elements of the 
deformation potential. In fact, the relative strength of the latter is strictly related to the 
bond stretching associated to the ionic displacements. In the present work, using the scheme 
introduced in ref. PU| for Ceo, we address the case of C70. 



II. ELECTRON-PHONON COUPLING 



Our calculation starts with the determination of the electronic states of fullerene-Cyo, 
and the corresponding electron density, within DFT-LDA. We include the exchange and 



correlation effects according to the parametrization of Perdew and Zunger [^] of the Monte- 
Carlo results of Ceperley and Alder No GGA corrections are considered, since LDA is 
known to yield very accurate results for the equilibrium geometry and charge density in such 
covalently bonded carbon sistems . The role of carbon atoms is taken into account using 
ab-initio norm-conserving pseudopotentials 0]. A spherical basis, made up with states 



which are solution of a central potential which mimics the L=0 component of the LDA local 
potential, is employed. Spherical wavefunctions with n and / respectively up to 15 and 20 
are included in the basis (this gives ~ 4000 basis vectors). The Kohn-Sham equations are 
solved in matrix form on this basis, as described in more detail in Ref. In that work 
the maximum values of n and / were slightly different, but the results are essentially the 
same. The HOMO-LUMO gap is 1.87 eV, in agreement with other LDA calculations [IB . 
As in Cgo, errors due to the neglect of both self-energy and excitonic effects nearly cancel 
each other, leading to a theoretical value close to the experimentally measured optical gap 
(1.6±0.2 eV). The Kohn-Sham levels up to about -10 eV are in overall agreement with the 
UPS spectra. The symmetries of the HOMO and LUMO levels {E'{) are the same as those 
obtained in 
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The different results of Ref. |26| are probably an artifact due to the limited 
basis set used there. 

The energies and eigenvectors of the vibrational normal modes have been taken from a 
BCM calculation Then, following pO|, we have parametrized the ab-initio electronic 

density in terms of sp^"*"^ hybrid orbitals obtained as linear combinations of the four s and 
p valence orbitals of each carbon atom. Three of the hybrid orbitals are directed along the 
bonds connecting the atom with its three nearest neighbors, while the fourth is determined by 
orthogonality conditions and takes care of the additional 7r-bonding present in the fullerenic 
cages. The radial wave functions of the s and p orbitals are taken as standard Slater wave 



functions {Rg 



(-r/o-a)^ 



where the parameters ai and (T2 have 



been adjusted in order to obtain the best fit to the C70 charge density computed in LDA. 
This leads to ai = 0.66 A and a2 = 0.36 A (to be compared with similar values obtained in 
the case of Cgo [13). For comparison, the atomic values are ai = 0.65 A and (72 = 0.17 A. 

After having determined the hybrid orbitals, we can write the contribution to the total 
electronic density arising from a single atom, and carry out a multipole expansion of it 
around the center of the cluster. Adding the contributions of the 70 atoms one obtains the 
total density. The deformation potential associated with a normal mode a is then 
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70 

E 

N=l 



QW ■ ^nVtot (r, {R}) 



(1) 



{R}={RO} 



where Qtv'' is the displacement of the iV-th ion, Vtot is the total LDA electronic potential, 
and {R} represents the whole set of ionic coordinates whose equilibrium values are {R^} 
(Qn = Rn — R-n)- The pseudopotential term of Vtot displays an explicit dependence on the 
ionic positions, therefore the calculation of its gradient presents no difficulties. On the other 
hand, the Hartree and exchange- correlation (XC) terms depend on {R} implicitly through 
the electronic density g. We thus write 



N 



Vat Vns,rtree,Xc[0{r,{R})] 



EdVua.rtTcc,Xc[Q] ^ 
A 



Vn Q{r,{R})- 



(2) 



One can easily calculate dV/dg, while the gradient of the total density is determined by 
assuming, within the present hybrid orbital model, that upon moving the ions according 
to the phonon eigenvectors the hybrid orbitals change their direction but not their shape. 
This makes the calculation quite simple and transparent, since it is possible to visualize the 
effect of a normal ionic displacement on the electronic orbitals. Finally, we can evaluate 
the matrix elements of the deformation potential, [V^'e/]jj = {i\V}^j\j), between Kohn-Sham 
wavef unctions. 

If we are interested in the behavior of the conductivity in alkali-doped fullerites, or in 
the photoemission spectrum of C70, we should start from phonon energies, and deformation 
potentials, evaluated by using the electronic density of a negative charged ion. On the other 
hand, because the density of the 280 valence electrons of C70 is not appreciably altered by 
adding a few more electrons, one expects the matrix elements of the deformation potentials 
associated with C70 and with C^q (1 < n < 4), to be rather similar. 

We evaluate, for each phonon a, the reduced matrix elements [V^e'/jij (corresponding to 
a given electronic state i) which are denoted by Qa, as well as the partial electron-phonon 
coupling constants Aq defined as 



XJN{0) = Cgl/hu, 



(3) 



In the above expression, A^(0) is the density of levels at the Fermi energy, uja are the phonon 
frequencies, and C is (ia/2, being the degeneracy of the phonon state |2^. The results 
corresponding to the LUMO {E'D and the A'^ and E'2 phonons, are collected in Table I. We 
only report those phonons yielding a non-negligible partial coupling, i.e., \a/N{fS) > 0.005 
meV. For these phonons, we also show the value of the full width at half maximum (FWHM) 
of the Raman peaks (7^), due to the decay to electron-hole pairs, widths which are connected 
to the partial couplings by the relation M 



7o 



(4) 



These widths show up as additional broadenings of the Raman peaks when the LUMO 
electronic state becomes occupied, e.g., in the case of K4C70. The differences between the 
widths of the Raman peaks of K4C70 and those of the peaks of bare C70 films have been 
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measured p9[, and are reported in the last column of Table I. However, a direct comparison 
with the theoretical results displayed in the fifth column of Table I is not straightforward 
for several reasons. First of all, it must be observed that subtracting the phonon widths 
measured in pure C70 from that measured in K4C70 leads, in some cases, to unrealistic 
negative values. Then, one should be aware of the fact that the frequencies and symmetry 
assignements of the experimental peaks are affected by some uncertainty. In any case, we 
followed the assignements given in refs. and []T^, which are different from those of ref. 
29| , since the former refer to low-temperature, higher resolution spectra, and rely also on 
calculations of the phonon peaks strength [^. Finally, the value of N{0) is not known with 
a high accuracy: we assume a value of 12 eV^^ as suggested by Ref. pOf, but any value 



within the range 6-30 eV~ seems to be equally possible fg]. The total electron-phonon 
coupling constant A (sum of the partial coupling constants) computed making use of the 
results reported in the fifth column of Table I, and of the empirical values of 7„ (seventh 
column of Table I), are quite similar and of the order of 0.1. Keeping in mind the above- 
mentioned caveats, we may hence speak of an overall agreement between the results of the 
present calculation and the experimental findings. If A^(0) varies between 6 and 30 eV~^, 
the theoretical total coupling constant varies between 0.05 and 0.25. These values must be 
compared with a value essentially equal to 1 obtained in the case of Cqo (cf., e.g., and 
pO| , as well as references therein). 

In a system where superconductivity is associated with the electron-phonon coupling, 
the transition temperature Tc can be calculated making use of McMillan's solution of the 



Eliashberg equation |3^ 



1.04(1 + A) 



A-/i*(l + 0.62A) 



(5) 



where cuin is a typical (logarithmic averaged) phonon frequency, and /x* the Coulomb pseu- 
dopotential ||33|. While the above expression leads to ~ 10-15 K for Cgo, making use of 
typical values fi* =0.2-0.3, it predicts a vanishing value for C70, even allowing for uncertain- 
ties of the order of 50 to 100 % in the estimated value of A. 

This result testifies to the fact that the transport properties in fuUerites are not simply 
dictated by the general features common to all these systems. Going from Cgo to C70, 
deformation plays a very important role with respect to the electron-phonon interaction, 
by breaking the symmetries of electronic and vibrational levels, and by weakening their 
coupling. 

In fact, in order to better understand the origin of the large reduction of the coupling 
strength in going from Cqq to C70, we have analyzed, for both clusters, the phonon eigenvec- 
tors of the low-frequency modes which couple to the electronic LUMO state, and compared 
them with the spatial localization of the latter. 



In the low-frequency region, as illustrated in ref. [|T^], it is possible to recognize, with no 
ambiguity, the relation between a group of slightly splitted vibrational modes of C70 and the 
corresponding quintuplet of an Hg mode of Ceo- This can be done by classifying the modes 
of both clusters according to the largest common symmetry subgroup, Cs^. We focus on the 
Hg{2) mode of Ceo at 408 cm~^, which is associated with one of the largest X^. Following 
ref. |]18[, we can trace the five modes associated with the Hg{2) vibration of Ceo to the A^,E^ 
and E2 modes of C70 at 425, 415 and 442 cm~^, respectively. The mode does not couple 
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with the LUMO of C70, hence about 2/5 of the couphng strength is lost. The remaining 
modes are associated with a higher Tiu with respect to the case of Ceo, further reducing the 
couphng. Finally, a large reduction effect comes from the different spatial localization of the 
LUMO wave functions in Cqo and C70. In particular, we have explicitly verified that the 
shape of the molecular deformation associated with the above mentioned normal modes is 
quite similar in the two clusters ||3^. However, in the case of Ceo the LUMO wave function 
is concentrated in the neighborhood of the bonds which undergo the largest stretching while 
in C70, the LUMO wave function is distributed on a larger region. We show in Fig. 1 the 
stereographic projection of Ceo LUMO wave function on a spherical surface with radius 3.4 
A, while in Fig. 2 a similar projection - but on an ellipsoidal surface whose radii are 3.55 A 
and 4.20 A - is presented for the case of the LUMO wave function of C70. In both figures, 
the thick continous lines mark the bonds which vary by more than 6% under the phonon 
displacement fields. 



III. PHOTOEMISSION SPECTRA 

The weakness of the electron-phonon matrix elements in the case of C70 reflects itself also 
in the shape of the photoemission spectrum of Cfo, as is shown below. In the photoemission 
experiment, we are dealing with the process C^^+hv C70 +e~, where hv is the incident 
photon and e~ the emitted electron. 

In order to define the ground state of Cfg, we start from the full electron-phonon Hamil- 
tonian, 

H = Y.e4c^ + Y.^uJJ:iv^ + Y.[vS]^3<^H^i + ^o.), (6) 

where the first term contains the Kohn-Sham energies e^, the second term is the free phonon 
term and the last term is the usual bilinear electron-phonon coupling term. The 280 valence 
electrons of C70 are assumed to be frozen (we call ipc their wavefunction), while the extra 
electron of Cfg is assumed to occupy one of the two (degenerate) lowest unoccupied states 
of C70 whose wavefunctions are labeled by (i=l,2). The basis used for solving Eq. (H) 
is built up with these one electron states coupled to 1, 2 or 3 different phonons (let us call 
0Q, a phonon wavefunction). The matrix thus obtained, which turns out to be of the order 
of 5000 x5000, is diagonalized by using the Lanczos method. We denote by £"0 its lowest 
eigenvalue and by \l/o the corresponding wavefunction. This should then describe the ground 
state of the charged molecule Cfo, 

'^0 = ^O'iAi'C+Yl ^»(") 0a + Yl 0a 0a' + ••• (7) 

Assuming that the emitted electron does not interact with the system left behind (sudden 
approximation ||3^), and that its wavefunction is described by a plane wave, the transition 
probability for the emission process, from the initial state \I'o to a final state of energy Ef in 
which the neutral C70 has 1, 2 or 3 phonons a excited (we denote by the corresponding 
total vibrational wavefunction), is given by the Fermi golden rule, 

W^~El(e"'''>c ^/iKxtl^o)!' 6{hu + Eo-Ef-e), (8) 
/ 
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where Kxt is the (dipole) external field, Ef = J2a (the energy of the 280 frozen electrons 
does not appear neither in the initial nor in the final state), and e is the electron energy. If we 
substitute Eq. (|^) into Eq. (|), and assume that the dominant terms are those corresponding 
to the LUMO (i.e., to the state ipi), we can factor out the dipole matrix element, and we 
are left with 

~ \ai\^ 6{hu + Eq - e) +J2 l^i(")r ^(hi" + Eq - - e) 

a 

+ ^ |ci(q;, a') Sihu + Eq — htUa — hua' — e). (9) 

We have plotted the results for Wq in Fig. 3. The scale has been set in such a way that 
hu + Eq = and the numbers in the horizontal axis correspond to e. While the shoulder 
at 300 cm~^ and the tail extending from 400 to 800 cm~^ indicate the presence of satellite 
peaks containing 1 and 2 phonons (as emphasized by the thin curves), the rather smooth 
behavior of the cross section and the strength of the main peak with respect to the satellites, 
as compared to the corresponding quantities in the case of Ceo, testify to the fact that the 
electron-phonon coupling in Cfg is much weaker than in Cgg (see inset of Fig. 3). While the 
overall behavior of the photoemission cross section predicted by theory is confirmed by the 
experimental data, lack of resolution does not allow for a detailed comparison. 

IV. CONCLUSIONS 

The calculated electron-phonon coupling in C70 is found to be about one order of magni- 
tude weaker than in Cgo- In particular, we estimate the total dimensionless electron-phonon 
coupling constant A to be about 0.1. The decrease with respect to Ceo is understood on 
the basis of the symmetry reduction, and the different shape and spatial localization of 
the LUMO wavefunction in C70. The calculated A value is consistent with the lack of the 
superconducting phase transition reported in the literature at temperatures as low as 1 K. 
This value agrees also with the overall features of two relevant experimental measurements: 
the broadening of Raman peaks in K4C70 and the photoemission cross section in Cfo- In 
the latter case, experiments with higher resolution would be required to test in detail the 
theoretical predictions. 
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TABLES 



TABLE L Reduced electron-phonon matrix element Qq, and partial coupling constants Xa/N{0) 
for the LUMO electron state and selected phonons of Cyo- The first two columns report the calcu- 
lated phonon frequencies and symmetries obtained within a bond-charge model calculation [^]. In 
the fifth column we also show the calculated widths of the phonon peaks (Eq. (^)). For comparison, 
the last two columns give the experimental frequencies and the mesured widening 7 of the Raman 
peaks in going from pure C70 to K4C70 films. 







Theory 






Experiment 
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Qa [meV] 


A«/iV(0) [meV] 
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1.08 


567 




3.5 
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3.48 
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25.7 
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6.79 
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30.4 
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FIGURES 



Fig. 1: Stereographic projection on a plane of the values taken by the squared Cgo 
LUMO wavefunction on a spherical surface of radius 3.4 A. The thick lines mark the bonds 
that vary more than 6% by effect of the displacement field associated with the Hg phonon 
at 408 cm~^ 



Fig. 2: Stereographic projection on a plane of the values taken by the squared C70 
LUMO wavefunction on an ellipsoidal surface of radii 3.55 A and 4.20 A. The thick lines 
mark the bonds that vary more than 6% by effect of the displacement field associated with 
the A'l phonon at 425 cm~^. 



Fig. 3: Calculated photoemission spectrum for Cfo (solid line). In the inset, the photoe- 
mission spectrum for Cgg from reference is shown. The contributions coming from zero 
phonons (long-dashed line), one phonon (dash-dotted line) and two phonons (dashed line) 
are also plotted separately. See the text for a discussion. A broadening with Gaussians of 
40 meV width has been used in constructing the photoemission spectra. 
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